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, ABSTRACT 

o . 

Aims. We calculate energy release associated with a first order phase transition at the center of a rotating neutron star. This 
. quantity is equal to the difference in mass-energies between the initial normal phase configuration and the final configuration 
containing a superdense matter core, with total baryon number and angular momentum kept constant. 

t Methods. The calculations of the energy release are based on precise numerical 2-D calculations, in which both the polytropic 

^ . equations of state (EOS) as well as realistic EOS of the normal phase are used. Presented results are obtained for a broad 
QO range of metastability of initial configuration and size of the new superdense phase core in the final configuration. When 
QQ the equatorial radius of the dense core of the superdense phase is much smaller than the stellar equatorial radius analytical 
r- H • expressions for the energy release are obtained. 

' Results. For a fixed "overpressure", SP, defined as the relative excess of central pressure of collapsing metastable star over the 
, pressure of equilibrium first-order phase transition, the energy release AE remarkably does not depend on the stellar angular 
momentum and coincides with that for nonrotating stars with the same SP. The energy release is proportional to (<5P) 2,5 for 
small SP, when sufficiently precise brute force 2-D numerical calculations are out of question. At higher SP, results of 1-D calcu- 
lations of AE(SP) for non-rotating stars are shown to reproduce, with very high precision, the exact 2-D results for rotating stars. 
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^ | 1. Introduction Pq. Consequently, a metastable core could form during 

•rH . _ neutron-star evolution in which central pressure increases, 

\S , One of the intriguing predictions of some theories of , , m, , 

/\ o r^ due to accretion or spm-down. 1 hen, nucleation of the cx- 

$H ' dense matter in neutron-star cores is a phase transition ,. , , ., s , . c ,. c c 

3 , . , , x °tic (higher density) phase implies formation ol a core of 

. >r into an exotic (i.e., not observed m laboratory) state. ,. , • • , , , , 

. t • • t i i • f • the exotic phase and is accompanied by a core-quake and 

Iheoretical predictions include boson condensation of pi- , . , . , , , , , , , 

, energy release. A theory which enables one to calculate 

ons and kaons, and dcconnncmcnt of quarks (sec, e.g., , , , i- j i c j. j 

, „ the changes m stellar parameters implied by a first-order 

Glendenning 2000, Weber 1999). , , 6 ... . p ... l \ { , , 

° ' phase transition m a non-rotating neutron star was devel- 

As far as the structure and dynamics of neutron stars Qped by Haensd et &L ig86 and zdunik ^ ^ lg8? (an 

is concerned, most important are the consequences of the e&rlier Newtonian theory was p reS ented by Schaeffer et al. 

first-order phase transitions accompanied by discontinu- 19g3) The energy released ^ & corequake was shown to 

ities in the thermodynamic potential densities. In the sim- depend strongly Qn ^ gize of thg denge ph&se corej thc 

plest case, one considers states consisting of one pure le&ding term bf ,. proportional to the fifth power of this 

phase. Because of high degeneracy of matter constituents, 

„ ,. i i T i i • core radius, 

effects of temperature are neglected, in thermodynamic 

. , ... , „ , In the present paper we calculate the energy release 

equilibrium, phase transition occurs then at a well denned . . . °^ 

, . . , n , ., . , , , due to a phase transition m a rotatinq neutron star. Our 

pressure Pq, and is accompanied by a density jump at the . f . . a . 

, . , r theory is based on the 2-D simulations and is much closer 

phase interface. , . , 

. „ , . ... .. . , ... to reality than the 1-D theory developed m (Haensel et 

A first-order phase transition allows for a metastabil- J t w 



ity of the pure "normal" (lower density) phase at P > 



al. 1986, Zdunik et al. 1987). In real world, evolutionary 
processes which lead to the increase of the central density 
Send offprint requests to: J.L. Zdunik in neutron star (accretion, slowing-down), as well as the 
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collapse itself, all take place in a rotating star. Of course, 
the 2-D calculations are incomparably more difficult than 
the 1-D ones. However, as we show in the present paper, 
when suitably parametrized, the energy released during a 
corequake essentially depends only on the excess of the 
central pressure of the metastable configuration over P , 
being to a good approximation independent of the angular 
momentum of collapsing star. 

The paper is organized in the following way. In Sect. El 
we introduce notations and describe general properties of 
the first-order phase transitions in stellar core with partic- 
ular emphasis on the metastability and instability of neu- 
tron star cores. Analytic considerations, concerning the 
response of a star to a first-order phase transition at its 
center, and in particular, the calculation , within the linear 
response approximation, of the energy release associated 
with such a transition, are presented in Sect. Analytic 
models of the EOSs with first-order phase transitions, al- 
lowing for very precise 2-D calculations, are considered in 
Sect. 13.11 where we derive generic properties of the en- 
ergy release due to a first order phase transition at the 
center of a rotating star. In Sect. 13.21 we present our re- 
sults obtained for a realistic EOS of normal phase, and we 
confirm remarkable properties of the energy-overpressure 
relation, obtained in the previous section. In Sect. [S] we 
present practical formulae suitable for the calculation of 
energy release associated with a first order phase transi- 
tion at the center of a rotating neutron star. Finally, Sect. 
contains discussion of our results and several examples 
of application of our formula for the energy release. 

2. EOS with a first order phase transition 

Let us consider a general case of a first-order phase tran- 
sition between the N (normal) and S (superdense) phases 
of dense matter. At densities under consideration, all con- 
stituents of matter are strongly degenerate, and temper- 
ature dependence of pressure and energy density can be 
neglected. At a given baryon density, rib, energy density 
of the N-phase of matter (including rest energies of parti- 
cles which are matter constituents) is fpj^b) and pressure 
•fk(^b)- The baryon chemical potential = enthalpy per 
baryon in the N phase, is /zn = (Pn + ^rO/ftb- Similarly, 
one can calculate thermodynamic quantities for the S- 
phase. 

A proper thermodynamic variable, continuous and 
monotonous in the stellar interior, is the local pressure 
P. Equilibrium state of the matter at a given P is real- 
ized at minimum of enthalpy per baryon. For P < Pq , this 
minimum is realized by the N phase, and for P > Po - by 
the S one. The value of Po is obtained from the crossing 
condition //n(P) = Ms(-P), which yields also the values of 
the matter densities, p^ and ps, and the corresponding 
baryon densities, tin and ns, at the N-S phase coexistence 
interface. These parameters are obtained assuming ther- 
modynamic equilibrium. A schematic plot of EOS of mat- 
ter with a first order phase transition N-S, in the vicinity 
of the phase transition point, is plotted in Fig.^ 
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Fig. 1. A schematic representation, in the p — P plane, of 
an EOS with a first order phase transition. Solid segments: 
stable N and S phase (in thermodynamic equilibrium). 
Dash-dot segment: metastable N phase. The S phase nu- 
cleates at Pnucl, which depends on the temperature and 
the compression rate. At P cr it nucleation of the S phase is 
instantaneous, because the energy barrier, separating the 
N phase from the S phase, vanishes. 

Solid segment of the N-phase curve, in Fig. ^ corre- 
sponds to stable N-phase state. For pressure above Po, the 
N phase becomes metastable with respect to the conver- 
sion into the S phase. The S phase can appear through 
the nucleation process - a spontaneous formation of the 
S-phase droplets. However, an energy barrier resulting 
from the surface tension at the N-S interface delays the 
nucleation for a time identified with a lifetime of the 
metastable state r nuc i. The value of r nuc i decreases sharply 
with P > Po, and drops to zero at some P cr it, where the 
energy barrier separating the S-state from the N-state van- 
ishes. For P > P cr ;t the N phase is simply unstable and 
converts with no delay into the S phase. 

Consider a neutron star built of matter in the N phase. 
Its central pressure P c increases during spin-down or ac- 
cretion. A quasistatic compression of the N-phase core 
moves matter into a metastable state with P > Pq. A 
metastable core of the N phase is bounded by a sur- 
face with P = Pq. As soon as the central compression 
timescale metastable core r comp = P c /P c becomes equal 
Tnuci, droplets of the S phase appear at the star center. 
This happens at central pressure P = P nU ci < Pcrit- The 
S-phase droplets introduce the pressure deficit, destabilize 
the metastable core and consequently - the whole star. The 
S-phase core grows, and the process of the N — >S phase 
transition proceeds until the final hydrostatic equilibrium 
state is reached. In the final state, a central core of the 
S-phase is bounded by an N-S coexistence surface, of con- 
stant pressure Po, on which matter density suffers a jump 
from p N on the outer N-phase side to p s on the inner S- 
phase side. 

In the present paper we restrict ourselves to the case of 
the density jump satisfying ps < |(p N +Po/c 2 ). Therefore, 
stellar configurations with a core of the S phase of arbitrar- 



J. L. Zdunik et al.: neutron-star core-quakes 



3 



ily small radius are stable (Seidov 1971, see also Kaempfer 
1981 and Zdunik et al. 1987). The case of a strong first 
order phase transition with pg > |(p N + Po/c 2 ), when 
configurations with a small S phase core are unstable and 
collapse into those with a large S phase core, will be pre- 
sented in a separate paper. 

A metastable stellar state of the N phase can be de- 
scribed by central overcompression - fractional excess of 
density relative to p N , 

S p = Ei—P^ . (i) 

Equivalently, we define a dimensionless central overpres- 
sure, 



Po 

Let us now consider the timescales and their inter- 
play. First, there is a microscopic nucleation timescale 
Tquci(^-P)- The time dependence SP(t) results from a global 
evolution of neutron star and depends on the astrophysi- 
cal scenario of central compression; it will be characterized 

by Tbomp = Pq/Pq = l/SP. The time needed to compress 
matter in the center of the star from Po to P c is, in lin- 
ear approximation, equal to r comp • SP. During this time 
the system remains in a metastable state. In contrast, the 
functional dependence T nuc \(SP) is given by the local dense 
matter microphysics and is a very sensitive function of 
P, being very large for P close to Po and dropping down 
abruptly above some pressure (see for example Iida & Sato 
1997). As a result, the condition for the time t = i qua ke 
at which the metastable stellar configuration collapses can 
be estimated by solving the equation: 

Tnucl [SP(t)] = Q ■ T comp • SP . (3) 

where a small dimensionless coefficient q reflects the very 
steep character of the function T nuc i(<SP) and is of the order 
10~ 2 — 10 -3 (details will be presented in our forthcoming 
paper). 

3. Calculation of the energy release 

We restrict ourselves to axially symmetric, rigidly rotating 
neutron stars in hydrostatic equilibrium. In what follows 
by "radius" we mean a circumferential radius in an equa- 
torial plane. 

We assume that at a central pressure P c = P nU ci the 
nucleation of the S phase in an overcompressed core, of ra- 
dius r N , of configuration C, initiates the phase transition 
and formation of S-phase core of radius r s in a new con- 
figuration C* , as presented on Fig. [3 Transition to an S 
phase at the core boundary, occurring at r s , is associated 
with a density jump characterized by A p = p s /p N . Because 
we are interested in relatively small S-phase cores, we will 
approximate its EOS by a polytrope with the exponent 
7 S , equal to the adiabatic index of the S phase at p = p s . 
Of course, in reality the adiabatic index of the S phase 



does depend on the density. However, for small cores the 
polytropic EOS is an excellent approximation, and in the 
limit when only leading terms in the r s -expansion are kept 
this approximation becomes exact (Bejger et al. 2005). 




Fig. 2. Transition from a one-phase configuration C with 
a meta-stable core of radius r N to a two-phase configura- 
tion C* with S-phase core with a radius r s . By "radius" 
we mean an equatorial circumferential radius. These two 
configurations have the same baryon number A = A* and 
total angular momentum J = J* . 

Having a pair of EOSs, one with and one without a 
softening by phase transition, the next step is to com- 
pare the hydrostatic equilibria of neutron stars corre- 
sponding to each of these EOSs. The models which we 
calculate are rigidly rotating, axisymmetric solutions of 
Einstein's equations. The numerical computations have 
been performed by means of a code built on the Lorene 
library (http://www.lorene.obspm.fr), with the accu- 
racy of 10~ 6 or better measured using general relativis- 
tic virial theorems. The neutron-star models can be la- 
beled by the central density p c and rotational frequency 
/ = n/2n. These parameters are natural from the point of 
view of numerical calculations (p c and / are input param- 
eters in the numerical code). But we can imagine another 
parametrization, more useful for other purposes. For ex- 
ample, to study the stability of rotating stars, the better 
choice is central density, p c , and total angular momentum 
of the star, J. 

Below p N , the EOS for two cases (with and without 
phase transition) is the same. For the problem of non- 
rotating stars, considered two decades ago papers (Haensel 
et al. 1986, Zdunik et al. 1987), the configuration with 
pc = P N was denoted by Co, and treated as a "reference 
configuration". Configurations Co, C, and C* are depicted 
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in the A — P c plane in Fig. |3| The radius, gravitational 
mass, and the total baryon number of Co will be denoted 
by Rq, Mo, and Aq. For rotating stars we do not have one 
" reference configuration" but a set of " reference configu- 
rations" which depend on rotation rate {Co(/)} or total 
angular momentum {Co (J)} of initial metastable config- 
uration. As it has been discussed in Sect. 2, there exists 
also a set of configurations resulting from the existence 
of metastability of the matter in N phase. The maximum 
central pressure which can be reached in the N-phase star 
is defined by the value P c — P nU ci- For P c > P nU ci, the 
phase transition in the center of the star takes place on a 
timescale much shorter than time of the stellar evolution 
(accretion or spin down). Thus the "critical line" defined 
by the condition P c — P nuc i corresponds to the configura- 
tions for which nucleation in the center triggers a collapse 
of the whole star and a corequake. ^From the astrophysical 
point of view, the phase transition in the center followed 
by a corequake takes place at a point of a "critical line". 
Consequently, "normalization" of results with respect to 
parameters of "reference configuration" is more compli- 
cated than in the case of non-rotating stars. There are 
many possibilities to define "reference configuration" for 
a star with metastable core (P c > P N ), rotating with fre- 
quency / and possessing a total angular momentum J and 
total baryon number A. Such a "reference configuration" 
should be defined as the configuration with the central 
pressure equal to Pq, as shown in Fig. [3] However, we can 
choose either a reference configuration star with the same 
J or the same / as that of collapsing metastable configura- 
tion C. In fact, to reach the metastable configuration with 
P c > Pq, we have to cross the "reference line" at the point 
defined by the evolutionary process leading to the insta- 
bility. Properties of "critical line" and "reference line" in 
the space of rotating configurations of slowing down or 
accreting neutron stars will be analyzed in more detail in 
our forthcoming paper. 

We additionally assume, that transition of the star 
from a one-phase configuration to the configuration with 
a small dense core built of S phase, takes place at fixed 
baryon number A (no matter ejection) and fixed total an- 
gular momentum of the star J (radiation loss of J ne- 
glected). 

The energy release during transition C{A 1 /) — > 
C*(A, /*) is calculated from the change of the stellar mass- 
energy during this process, 

AE=[M{C)-M{C*)] A<J (? (4) 

3.1. Energy release for polytropic EOSs 

In the present section we intend to simplify the analysis, 
removing problems connected with numerical precision. 
To this aim, we will use the polytropic EOSs for the N 
and S phases. The polytropic EOSs not only guarantee 
high precision of numerical calculation, but also open a 
possibility of the exploration of wide region of the param- 
eter space. Description of the polytropic EOSs and their 
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Fig. 3. Total baryon number A of hydrostatic stellar con- 
figuration, versus central pressure P c , at fixed stellar an- 
gular momentum J, for the EOS, depicted in Fig.^ Solid 
line denotes stable , dash-dot line - the states which a 
meta-stable with respect to the N — >S transition. For a 
central pressure P nU ci the S-phase nucleates in the super- 
compressed core of configuration C, and this results in a 
transition C — > C* into a stable configuration with a S- 
phase core and central pressure P* . Both configurations C 
and C* have the same baryon number A. 

application to relativistic stars with phase transitions was 
presented in detail in our previous publications in this se- 
ries (Bejger et al. 2005, Zdunik et al. 2006). 

In Fig. 0] we presented the energy release as a function 
of r s = r s /R{C), for several values of the angular mo- 
mentum of the metastable configurations C(J). The ref- 
erence configuration is C - the initial configuration of the 
corequake with metastable core in the center and central 
pressure P c = P nU ci- As we see in Fig. 0J the energy re- 
lease corresponding to a given value of r s depends rather 
strongly on the rotation rate (here presented for the fixed 
values of total angular momentum). 

In Fig. we presented the energy release as a function 
of the overpressure of the metastable N phase in the cen- 
ter of the metastable star C(A, J), SP = P nU ci/Po — 1, for 
several values of J. As we already stressed, the value of 
Pmci (or SP) can be determined from microscopic consid- 
erations, combined with physical conditions prevailing at 
the star center as well as with their time evolution rate. 
Having determined P nU ci, we can determine the energy re- 
lease, AE, due to the corequake C(A,J) — ► C*(A,J), 
where the metastable one-phase configuration, and the fi- 
nal two-phase configuration, have the same values of the 
baryon number A and total angular momentum J, which 
are conserved during the transition. 

As we see in Fig. [31 the energy release in a 
(mini)collapse of rotating star is independent of rotation 
rate of collapsing configuration, and depends exclusively 
on the degree of metastability of the N phase at the stellar 
center (departure of matter from chemical equilibrium), 
measured by the overpressure SP. In particular, to calcu- 
late the energy release associated with a corequake of a 
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Fig. 4. The energy release due to the corequake of ro- 
tating neutron star as a function of the dimensionless 
equatorial radius of the new phase core, r s . Different 
curves correspond to the different values of total angu- 
lar momentum of rotating star, fixed along each curve, 
J = (0,0.1, . . . , 0.9) x GM|/c, from the bottom curve 
to the top curve. 

rotating neutron star, it is sufficient to know the value of 
AE for a non-rotating star of the same central overpres- 
sure, which can be deduced from expressions derived in 
Zdunik et al. (1987)| 

It should be stressed that the configurations C(A, J) 
and C*(A,J), considered in this section, are really fast 
rotating ones, not so far from the Keplerian limit. This is 
visualized in Fig. 13 where we plotted the oblateness of the 
star and the kinetic to potential energy ratio. And still, 
in spite of fast rotation and large oblateness, the energy 
release is the same as in a non-rotating star of the same 
initial central overpressure. 

3.2. Energy release for realistic EOS 

In the present section we consider realistic EOS of the N 
phase. In order to explore how a realistic rotating neutron 
star will respond to the appearance of a core of S phase 
core, we used a recent SLy EOS of Douchin & Haensel 
(2001). The SLy EOS describes in unified way (i.e., start- 
ing from a single effective nuclear Hamiltonian) both the 
crust and the core of neutron star. In its original version, 
the SLy EOS assumes that neutron star core is composed 
of neutrons, protons, electrons and muons. We introduced 
a softening by a first order phase transition at tin, as 
shown in Fig. [7] 

We considered two EOSs of the S phase. Both were 
of polytropic form P(nb) = K s n^ s , with 7 S = 2.5 and 



Fig. 5. (Color online) The energy release due to the mini- 
collapse of rotating neutron star as a function of the over- 
pressure 5P of the N phase of the matter in the center 
of the star, for the polytropic EOS (7 = 2) with first 
order phase transition located at three different points 
(pressures Pq) for which the mass of the reference non- 
rotating configuration is equal to 0.7, 1.0, 1.2M Q . For 
each value of P two models of S phase EOS are con- 
sidered (one of examples of the EOS depicted in Fig. 0). 
They correspond to the S phase described by the poly- 
tropes with 7 S = 2.5 (upper curve) and 7 S = 3 (lower 
curve). The points of different color correspond to the dif- 
ferent values of total angular momentum of rotating star, 
J = (0, 0.1, . . . , 0.9) x GM%/c. For a given EOS, results 
for all rotating configurations can be very well approxi- 
mated by a single curve, independent of J. 

7 S = 3, respectively. We assume the values of n N and 
n s = A n n N . The values of n N and X n determine then Pq. 
Equality of baryon chemical potentials at P = Pq fixes the 
other constant of the S-phase EOS, which is the value of 
energy per baryon (including rest energy) at zero pressure. 
In this way, the values of n N and A„, together with 7 S , fully 
determine the EOS of the S phase. 

In Fig. 03 we show the energy release due to the 
C(A,J) — > C*(A,J) transition, versus overpressure, for 
two selected values of 7 S . The values obtained for different 
values of J are marked with different color and symbols, 
and can be distinguished on the electronic versions of this 
paper. However, as in the case of the polytropic models of 
the N phase, all color points lie along the same line. For a 
given overpressure SP, the energy release does not depend 
on J of collapsing metastable configuration, confirming in 
this way results obtained for the polytropic EOSs. This 
property is fulfilled astonishingly well for a broad range of 
of stellar angular momentum, J = (0.1, . . . , 0.7) x GMq/c. 
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Fig. 6. Top panel: the ratio of polar radial coordinate to 
the equatorial radial coordinate radius versus frequency 
of rotation. Bottom panel: the ratio of the kinetic energy, 
T and the absolute value of the potential energy, W ver- 
sus stellar angular momentum. Calculations performed for 
stellar configurations consisting of the N phase of dense 
matter, described by the polytropic EOS. Large dots cor- 
respond to the values of the total stellar angular momen- 
tum, J = (0.1, ... , 0.9) x GMq/c, which were used in Figs. 

HE] 

Our numerical results for different equations of state 
and models of a phase transitions presented in Figs [SJ |H1 
show that the energy release associated with a phase tran- 
sition in the core remarkably does not depend on the stel- 
lar rotation rate for a fixed overpressure at which the phase 
transition takes place. To estimate the 'accuracy' of this 
conclusion we can calculate the departures of the energy 
release for the rotating stars from the curve determined for 
nonrotating configurations. The relative differences are of 
the order of 1-2% for 5P > 0.1 (except for the smallest 
masses in Figs El 00 (M = O.7M and M = 1M© respec- 
tively) where they can be as large as 5%). For smaller 
SP > 0.1 the relative differences are larger because the 
absolute value of AE is very small and the differences are 
comparable to the numerical accuracy. 

4. Calculation of the energy release using linear 
response theory 

Below p N , the EOS for two cases (with and without phase 
transition) is the same. The ambiguity of definition of Co 
for rotating stars was discussed in Sect. 3. We can choose 
a reference configuration star with the same J or the same 
/ as that of collapsing metastable one. In fact, to reach 
the metastable configuration with P c > P N , we have to 
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Fig. 7. The EOSs with first order phase transitions, used 
in the present paper. Upper panel: the SLy EOS for the N 
phase, with a phase transition to the dense S phase. We 
assume n N = 0.54 fm -3 and A n = n s /n n = 1.25, which re- 
sults in A p = p a /p N = 1.28987. For the further description 
of the EOS, see the text. Lower panel: polytropic models 
of the EOS with a first order phase transition. We assume 
n N = 0.25 fm~ 3 and A„ = n s /n^ — 1.25, which results in 
A p = p s /p N = 1.2647. For the further description of these 
EOSs, see the text. 



cross the "reference line" at the point defined by the evo- 
lutionary process leading to the instability. 

We assume, that transition of the star from a one-phase 
configuration to the configuration with a small dense core 
built of S phase, takes place at fixed baryon number A 
(no matter ejection) and fixed total angular momentum 
of the star J (radiation loss of J neglected). Thus for a 
fixed baryon number A and fixed J, we will calculate the 
mass-energy difference between C and C*. 

We assume that the mass of the S-phase core is much 
smaller than the total stellar mass. Therefore, from now on 
we will be able to restrict ourselves to a linear response of 
neutron star to the appearance of the new dense S-phase 
core. The calculation is based on expressing the change in 
the density profile, due to the presence to a small core, 
as the combination of two independent solutions of the 
linearly perturbed equations of stellar structure (Haensel 
et al. 1986, Zdunik et al. 1987). The presence of a denser 
phase in the core changes the boundary condition at the 
phase transition pressure Pq and allows us to determine 
the numerical coefficients in the expression for the density 
profile change. The leading term in the perturbation of the 
boundary condition at the edge of the new phase results 
from the mass excess due to the higher density of the new 
S phase as compared to the N one in the supercompressed 
core. 
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Fig. 8. (Color online) The energy release due to the mini- 
collapse of rotating neutron star as a function of the 
metastability of the normal phase of the matter in the 
center of the star for the SLy EOS with first order phase 
transition located at two different points (pressures Po) 
for which the mass of the reference nonrotating configu- 
ration is equal to 1.0, 1.4M Q . For each value of Po two 
models of the dense core EOS are presented - polytropes 
with 7 = 2.5 (upper curve) and 7 = 3. The points of 
different color correspond to the different values of total 
angular momentum of rotating star. The results for all 
rotating configurations can be very well approximated by 
one curve. 

4.1. Non-rotating neutron stars 

Let us remind expressions obtained for static case consid- 
ered by Haensel et al. (1986), and Zdunik et al. (1987). 
For a non-rotating configuration of hydrostatic equilib- 
rium, the density jump p N — > p s leads to the lowest-order 
expression for the core-mass excess (with respect to the 
pure N-phase configuration), 



5 m 



core ., ~(/' : . 



0(r 5 



(5) 



The radius of the overcompressed N-phase core, r N , is 
connected with the overpressure which can be achieved by 
the N phase of matter, by 



5P=?lGx„p N (l+x n )(l + 3x„)rl . (6) 

The conservation of the stellar baryon number during the 
C — ► C* transition, implies a relation between r N and r s 
which neglecting terms r 4 and higher, reads 



(1 + 3xJr 2 N = (3 - 2A p + 3x N )r 2 B - (1 - A p ) a N r 3 s 
where a N is a numerical coefficient depending on Co 



(7) 



Using then Eqs. ®) and 0, one concludes that the 
leading term for the energy release expression oc con- 
while the next 
- I) 2 . The final 



tains a prefactor (A p — 1)(3 — 2A p + 3x h 
order term oc r\ 



contains a prefactor (A p 



expression for a normalized energy release reads 
M - M* 



AE = 



M 
+ (A P - I) 2 



(A p -l)(3-2Ap + 3z N )ai r s b + 



a 2 r s 



(8) 



where r s = r s /R . The coefficients ot\ and a-i are func- 
tional of the static reference configuration, Co (J = 0). 
Expressing now r s in terms of overpressure, we get the 
energy release in terms of the overpressure, 



AE- 



50 



(A P - 1) 



(3-2A p + 3;r N ) 1 - 



r/3i (SP) 



~D\2.5 



(9) 



where we use standard notation P50 = E/10 erg. The 
coefficient /3i depends only on the EOS of the N phase. 
More precisely, this coefficient is a functional of the non- 
rotating reference configuration Co and depends rather 
weakly on the mass Mq of the reference configuration 
Co for a rather broad range of masses, between 1M Q 
and 0.8M max . For EOSs considered in the present pa- 
per, numerical calculations give in a good approximation 
f3 1 ~ 0.016. 



4.2. Rotating neutron stars 

Numerical results of Sect. 3 clearly show that for for a 
given central overpressure of collapsing configuration, SP, 
the energy release does dot depend on J and coincides 
with that for non-rotating stars. Therefore, the formula 
© of Sect. 14. ll is valid also in the rotating case, provided 
<5P is sufficiently small. Our calculations show that this 
formula is quite precise for SP ?S 0.05 . The prefactor in 
front of [3\ , involving A p and x N , is identical with those 
obtained for the J = transition in (Haensel et al. 1986, 
Zdunik et al. 1987). The coefficient /3i depends only on 
the EOS of the N phase. More precisely, this coefficient is 
a functional of the nonrotating reference configuration Co- 

5. Numerical estimates of the energy release 

For practical application, and for small overpressures, 
SP < 0.1, it is convenient to summarize results obtained 
in Sects. 1^113. 21 in a formula 

AP 50 = a a (SP) 2 5 (10) 

The only dependence of coefficient a\ on the phase tran- 
sition parameter A p is via prefactors, 



(A P - 1) 



(3 - 2A P 



3x, 



(11) 



which allows for a rapid re-calculation of a\ when one 
changes the value of the density jump. 
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For the phase transition model considered in Sects. HOI 
the energy accompanying phase transition in a metastable 
star with central overpressure SP = 0.1 is about 2 x 
10 50 erg, and becomes one order of magnitude smaller for 
SP = 0.05 

6. Discussion and conclusions 

The most important result of the present paper is that the 
total energy release associated with a first order phase 
transition at the center of a rotating neutron star does 
depend only on the overpressure at the center of the 
metastable configuration and is independent of the star ro- 
tation rate. This result holds even for fast stellar rotation, 
when the star shape deviates significantly from spheric- 
ity, and J ~ 0.9GM|,/c, and for overpressures as high as 
(10-20)%. The property is of great practical importance, 
because it implies that the calculation of the energy re- 
lease for a given overpressure, requiring very high preci- 
sion to guarantee A = A*, can be reduced to the case of 
non-rotating spherically symmetric stars. 

To illustrate the practical importance of our result, let 
us consider the case of overpressure SP = 0.03, when the 
energy release is 2 x 10 49 erg. For a 1.4 Mq star this con- 
stitutes about 10~ 5 of the stellar mass-energy. Therefore, 
to get a meaningful result for the energy of this magni- 
tude, conservation of the total baryon number has to be 
satisfied at the level of one part in a million, which is 
next to impossible for a 2-D calculation with a realistic 
EOS. However, such a precision can be easily reached for 
spherical stars, provided the EOS is used in a thermody- 
namically consistent way (see, e.g., Haensel & Proszynski 
1982). 

Moreover, we have shown that for overpressures 
smaller than 5%, the energy release is proportional to 
power 2.5 of the overpressure, with coefficient weakly de- 
pendent on the mass of collapsing configuration, for stellar 
masses in the range from 1 M to about 0.8 M max . 

The energy release AE that we calculated is an ab- 
solute upper bound on the energies which can released 
as a result of a phase transition at the star center. The 
available channels may include, for example, stellar pul- 
sations, gravitational radiation, heating of stellar interior, 
and X-ray emission from neutron star surface. Moreover, 
the phase transition in rotating star implies shrinking of 
stellar radius, decrease of moment of inertia, and spin-up 
of rotation. These topics will be considered in our next 
paper. 
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